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ABSTRACT 

We derive the evolution equations for the electric and magnetic parts of the 
Weyl tensor for cold dust from both general relativity and Newtonian gravity. In 
a locally inertial frame at rest in the fluid frame, the Newtonian equations agree 
with those of general relativity. We give explicit expressions for the electric 
and magnetic parts of the Weyl tensor in the Newtonian limit. In general, the 
magnetic part does not vanish, implying that the Lagrangian evolution of the 
fluid is not purely local. 



Subject headings: cosmology: theory — large-scale structure of the universe 
gravitation 



1. Introduction 



Gravity is a long-ranged force. According to this Newtonian perspective, the motion of 
a mass element is affected by distant elements. In general relativity, however, the motion 
of a freely-falling body is determined by the local curvature of the spacetime manifold: in 
a locally inertial frame, the motion of nearby mass elements is governed by the Riemann 
curvature tensor. The Newtonian and relativistic viewpoints are made consistent by the 
fact that the Riemann tensor incorporates Newtonian gravitational tidal fields. 
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In the Newtonian approach, tidal fields are obtained from the gradient of the gravity 
vector, whose determination requires a sum over all mass elements. Thus, it came as 
a surprise to us that, in general relativity, the tidal field evaluated at the position of a 
freely-falling mass element might evolve according to a purely local equation — provided 
that certain stringent conditions are met. 

Matarrese, Pantano, & Saez (1993) and Croudace et al. (1994), following earlier work 
of Barnes & Rowlingson (1989), showed that if a quantity known in general relativity as the 
magnetic part of the Weyl tensor vanishes in the Newtonian limit, the Newtonian tidal field 
obeys a local Lagrangian evolution equation until trajectories intersect. In other words, the 
tidal field following a mass element changes with time in a way depending only on local 
fluid variables — the density and velocity gradient. Because these variables themselves 
obey Lagrangian equations (local aside from the tidal field) — it seemed that the one might 
be able to evolve the fluid variables and tidal field independently for all mass elements until 
trajectories intersect. In this situation, all nonlocal information is incorporated into the 
initial value of the tidal field. 

Bertschinger & Jain (1994) used this fact to study the nonlinear evolution of density 
perturbations in the expanding universe. They noted (as did Matarrese et al. 1994) that 
the magnetic part of the Weyl tensor did not necessarily vanish in the Newtonian limit. 
However, they assumed that it could be approximately neglected, and solved the coupled 
nonlinear fluid and tidal evolution equations for general irrotational motion starting from 
the growing mode of cosmic density perturbations. They found that nonlinear coupling 
of the fluid shear and tidal field would favor filamentary gravitational collapse as opposed 
to the sheetlike pancake collapse predicted on the basis of kinematical theory (Zel'dovich 
1970). This surprising conclusion rests on an important unchecked assumption. 

Until now there has been no Newtonian derivation of the tidal evolution equation. 
The equations of motion for the Riemann tensor in general relativity have, instead, been 
projected into the local fluid frame. This work was pioneered by Kundt & Trumper (1961), 
Hawking (1966), and Ellis (1971, 1973). It has resulted in a powerful covariant Lagrangian 
fluid description of matter and gravitational fields. This method has been applied to the 
evolution of cosmic density fluctuations by Hawking (1966), Ellis & Bruni (1989), Hwang 
& Vishniac (1990), and many later workers. In linear perturbation theory the magnetic 
part of the Weyl tensor vanishes for irrotational perturbations so that the Lagrangian fluid 
equations reduce to local equations. In cosmological perturbation theory the Lagrangian 
fluid approach, while elegant and free of gauge ambiguities, offers no compelling advantage 
compared with traditional Eulerian methods. 

Traditionally, two approaches have been used to study gravitational dynamics in 
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cosmology. The first is Eulerian: the fluid and gravitational variables (mass density, velocity, 
etc.) are defined on a grid and evolved according to partial differential equations. This 
method works well for a collisional fluid in which pressure forces prevent the intersection of 
trajectories, but it breaks down for cold dust (pressureless collisionless matter, e.g., cold 
dark matter) after mass elements intersect. The second method is that of Lagrangian 
trajectories, exemplified by N-body simulations: the mass is discretized into particles whose 
positions and velocities are integrated using the appropriate equations of motion. 

The Lagrangian fluid approach offers a third way: the density and velocity gradients, 
but not necessarily the positions and velocities, are integrated for individual mass elements. 
Clearly, this procedure is incomplete without computing the trajectories of mass elements. 
One can state the density of a given mass element (provided that it is not superposed 
on another element) but, without integrating the trajectories, one cannot say where it 
is. While the trajectories can be integrated here as in N-body simulations, this presents 
an extra complication. For some purposes, we would be happy to know the density and 
velocity gradient of every mass element even if their positions are unknown. 

This explains the attraction of local Lagrangian evolution. However, the applicability 
of the method remains unclear as long as questions remain about the magnetic part of the 
Weyl tensor. 

In this paper we address three questions: Does the magnetic part of the Weyl tensor 
vanish in the Newtonian limit? Can the Lagrangian evolution equations for the Weyl tensor 
be derived in the Newtonian limit? Are they local in general? In the following sections we 
show that the answers are, respectively, not necessarily, yes, and no. 



In this section we present the general relativistic derivation of the evolution equations. 
Although this material largely repeats the work of Ellis (1971, 1973), it is a necessary 
prelude to the Newtonian derivation. Besides defining the Weyl tensor and other quantities, 
we establish some notation that is used throughout this paper and we clarify the conditions 
required to obtain the Newtonian limit of general relativity. We adopt the conventions and 
notations of Misner, Thorne, & Wheeler (1973), including metric signature +2. 

The Weyl tensor is the traceless part of the Riemann curvature tensor (Misner et al. 



2. Evolution of the Weyl tensor in General Relativity 



1973): 
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where, for convenience in what follows, we have defined 

(2) 

The full Riemann tensor follows from the Weyl tensor and the Ricci tensor R^ v = R K ^ KU 
and its trace R = i? M „. Through the Einstein field equations, 

-R^ — -^Qfiu = , (3) 

the Ricci tensor gives the contribution to the spacetime curvature from local sources with 
energy-momentum tensor T^. The Weyl tensor gives the contribution due to nonlocal 
sources. Therefore, Newtonian tidal forces will be represented in the Weyl tensor. 

In the most common formulation of general relativity, the Einstein equations are 
regarded as field equations for the metric tensor components, since the Ricci tensor can be 
obtained from the metric tensor and its first and second derivatives. Because the Riemann 
tensor also follows from the metric and its derivatives, the Weyl tensor can be calculated 
from the solution to the field equations. 

However, we are interested in an alternative formulation of general relativity in 
which the Weyl tensor is treated as the fundamental geometrical quantity and the Ricci 
tensor follows algebraically from the Einstein equations for a given distribution of energy 
and momentum. This formulation requires equations of motion for the Weyl tensor 
independently of the Einstein equations. They follow from the contracted Bianchi identities 
(Kundt & Trumper 1961). In terms of the Weyl and Ricci tensors, these are 

^ K C^ UK \ = V[^R u ]x + - g\[ fl S/ u ]R K K . (4) 

Gradient symbols denote the covariant derivative with respect to g^ u . Square brackets 
around a pair of indices denote antisymmetrization, e.g., A^ = ^(A^ — A Ufl ). Symmetrized 
indices are surrounded by parentheses: A^ = ^(A^ u + A Ufl ). Substituting the Einstein 
equations (^]) into the Bianchi identities (Q) now provides field equations for the Weyl 
tensor: 

V K C^ KX = 87rGfa lfl T u]x + ~g x[ ^ v] T K l ^ . (5) 

Although the Weyl tensor has 256 components, only 10 of them are independent in 
four dimensions. It is convenient to incorporate these 10 components into two symmetric 
second rank tensors using the 4- velocity field u^(x) to split the Weyl tensor as follows: 

E^(u) = u*u x C^ uX , H^{u) = X - u K u x e aMll C a(3 u)x . (6) 
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These tensors, called the electric and magnetic parts of the Weyl tensor, respectively, fully 
determine the Weyl tensor for any non-null normalized m m (not just the 4- velocity of the 
matter). Indeed, the full Weyl tensor can be reconstructed from its electric and magnetic 
parts (Ellis 1971): 

C^ukX = g K \-y5 - t/iuap e K \ys) u a u f E ps (u) + (e^ Q/3 g K \ l5 + g^p e^x-ys) u a u 1 H ps {;u) . (7) 

Equation (|7|) is the exact inverse of equations (Q) provided g^u^v? = ±1 and the same 
is used in both equations. Ellis (1971) has a sign error in the first term of his version of 
equation (|7|) at the end of his section 4.3.3. Note that we have used the fully antisymmetric 
tensor e^ VK \ = (—g) 1 ^ 2 [pvnX], where g is the determinant of g^ v and [punX] is the completely 
antisymmetric Levi-Civita symbol defined by three conditions: (1) [0123] = +1, (2) [pvnX] 
changes sign if any two indices are exchanged, and (3) [pvnX] = if any two indices are 
equal. (Ellis uses the tensor r]^ lUK \ = —e^ K \. We have compensated for the sign change in 
defining H^.) 

The Weyl tensor is completely independent of the 4-velocity u^. Only the electric and 
magnetic parts depend on — in particular, they are orthogonal to E^ v u u = H fiu u u = 0. 
Thus, in the local rest frame defined by u^, only the spatial components of E^ u and H^ u 
are non- vanishing. Moreover, these tensors are trace-free: = H^^ = 0. Each therefore 
has 5 independent components, accounting for the 10 independent components of the Weyl 
tensor. 

Equation (||) yield equations of motion for the electric and magnetic parts of the Weyl 
tensor. The results depend, of course, on the energy-momentum tensor T^ v . We assume 
that the matter is a perfect fluid for which T^ u = (p + p)u fl u v + pg^ v with u M being the fluid 
4-velocity and p and p being the proper mass density and pressure, respectively, in the fluid 
frame. Using the same 4-velocity to split the Weyl tensor into its electric and magnetic 
parts, we obtain the following equations (Hawking 1966; Ellis 1971): 

(div-E) : P^P^V v E a(} + e^u u a ai H\ - 3H^ = ^ GP^V u p , (8) 

(H) : P^g^EL - P^e^UpV^s - 2u a apE^e u)a ^ + QH^ 
dX 

+P^(o aP H aP ) - 3H a ^a v) a + H a ^oo u) a = , (9) 
(div-if) : P» a P v pV v H af ' - e^ aP u v a ai EP p + 3E» v uj u = -8nG(p + p)^ , (10) 

(E) : p^ E^- + P^e^upV^Has + 2u a apH + 0£T 

dX 

+P^(a a ^E aP ) - 3E a ^a u) a + E a ^u v \ = —AtcG(p + p)a^ v .(11) 

We have adopted Ellis' names for these equations, inspired by their similarity to the 
Maxwell equations (Ellis 1971, 1973). 
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In writing equations fl8|)— (|TTD we have defined several new quantities. The derivatives 
of E a @ and H af3 are projected into the rest space of u v using the projection tensor 
Pfiu = g^u + UfjU v . The covariant derivative along the fluid worldline is D/dX = u u V u . The 
velocity gradient has been decomposed into the acceleration 4- vector a„, expansion scalar 
9, shear tensor cr^, and vorticity tensor u)^ or its dual uj^\ 



6 = V^m m , cr^ = a^ v ) , = = e^afsu ^ . (12) 



In the fluid rest frame, uJi is half the usual three-dimensional vorticity (V x v)^ Ellis (1971) 
defines and uj^ with the opposite sign. 

Equations @ and (|TT|) represent Lagrangian evolution equations for the Weyl tensor: 
the time derivatives are taken along the fluid worldlines. Equations (§) and (|II]) are 
constraint equations. However, only the divergences of E^ v and are constrained; the 
"curl" parts (the gradient terms in eqs. || and |11|) may be specified arbitrarily on some 



initial hypersurface. This is in contrast with Newtonian gravity, where the gravitational 
field is fully determined from the matter distribution through the static Poisson equation. 
It is similar to electromagnetism in that the gravitational field may contain a source-free 
part corresponding (in the weak-field limit) to gravitational radiation. 

Equations flSD-flTTl) are fully covariant tensor equations and can be applied in any 
coordinate system. Because of the key role played by the fluid 4- velocity, however, it is most 
convenient to evaluate them in the fluid rest frame. We shall assume that the matter is 
cold (pressureless) dust, in which case the 4-acceleration vanishes, so that the fluid rest 
frame is a locally inertial frame. Using locally flat coordinates in this frame, the equations 
of motion for the electric and magnetic parts of the Weyl tensor become 

Vj&i - e ijk a jl H\ - ZHijJ = y GV lP , (13) 

^ + V fc e kl {i E j)t + QH l3 + Sy a kl E kl - 3a k {i H j)k - u k {i H j)k = , (14) 

VjH^ + e ijk a jl E k l + SE^u? = -torGpuj , (15) 

dE- ■ 

—± - V k e kl {i H j}l + GE i:j + Sa a kl E kl - 3a k (i E j)k - tu k (i E j)k = —AnGp a tj . (16) 

Note that these equations assume w M = (1,0) at the point where they are being applied 
because we are using locally flat coordinates in the fluid rest frame, but the gradient of the 
3-velocity does not necessarily vanish. In section 4 we shall furnish a Newtonian derivation 
of these equations. 
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3. Weyl Tensor for a Perturbed Robertson- Walker Spacetime 

As we saw in the preceding section, the electric and magnetic parts of the Weyl tensor 
are only partly constrained by the matter distribution. They are fixed, however, when the 
metric is specified. In this section we shall obtain the Weyl tensor components and its 
electric and magnetic parts for a perturbed Robert son- Walker spacetime. 

We start with the following line element: 

ds 2 = a 2 (r) {-(1 + 2tp)dr 2 + 2w i drdx l + [(1 - 20) 7ii + dx'da?} , j ij h i:j = . (17) 

The perturbations ip, 0, i^i, and h^ vanish for a Robert son- Walker spacetime with expansion 
scale factor a(r) (r is conformal time) and 3-metric jij (the metric of a constant curvature 
space). We treat the perturbations as being small quantities. 

As written, our metric is completely general. To reduce the gauge freedom we impose 
the following transversality constraints on Wi and h^: 

ryVViWj = , y*V*/ly = , (18) 

where Vj is the spatial covariant derivative relative to 7^, whose inverse is 7 tJ ". It can be 
shown that, with these conditions, ip and are identical to the gauge- invariant scalar mode 
variables $a and — respectively, of Bardeen (1980). Moreover, u>j corresponds to the 
"vector mode" (gravitomagnetism) and h^ to the "tensor mode" (gravitational radiation). 
As discussed in detail by Bertschinger (1994), the gauge choice implied by equation (|i~8f) has 
several advantages over other choices such as the synchronous gauge. First, it is essentially 
unique — there is no residual gauge freedom associated with spatially inhomogeneous 
redefinitions of the coordinates. Second, by eliminating scalar mode contributions from Wi 
and hij, and vector mode contributions from h^, the physics of the modes is simplified. 
Finally, the perturbation variables are small provided that physical curvature perturbations 
are small, so that the coordinates (r, x l ) are nearly identical to locally flat coordinates 
scaled by a homogeneous expansion factor. These last two advantages facilitate relating 
calculations in this gauge to the Newtonian limit (i.e., locally flat spacetime with slow 
source speeds relative to the comoving frame). 

Before giving the Weyl tensor components, we first present the equations of motion for 
the metric perturbation variables obtained from the Einstein equations. These will provide 
intuition about the physics of the different types of perturbations (scalar, vector, tensor) 
and will facilitate comparison with the Newtonian limit. 

The scalar and vector metric perturbation fields obey the following equations 
(Bertschinger 1993) derived from the Einstein equations in the limit of distance scales small 
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compared with the curvature and Hubble distances and with nonrelativistic shear stresses 
(the latter condition implying if) = 0): 

V-g = -AirGa 2 (p-p) , Vxg + d T H = 1 V • H = , V x H = -lQirGa 2 f ± , (19) 

where 

£ = -VV> - <9 r w , H = V xw . (20) 

We use vector notation for three-dimensional vectors in the spatial hypersurfaces with 
metric 7y. Note that p(r) is the background density and f± = (pv)± is the transverse 
energy current (transversality implying V • fx = 0) evaluated in the comoving frame (so 
that v is the peculiar velocity, v = dx/dr to lowest order in the metric perturbations). 



If one neglects w, then equations (p!9|)-(p0|) reduce to Newtonian gravity in comoving 
coordinates. Because the source for the gravitomagnetic field H is smaller by 0(v/c) than 
the source for Newtonian gravity, H usually is unimportant. Note that equations (|T9| ) are 
nearly identical with the Maxwell equations; the important difference is the absence of 
longitudinal current and a displacement current (d T g ) term in the "Ampere" law for V x H . 
This difference implies that both g and H are essentially static fields without radiation. 

The absence of gravitational radiation from g and H is what one expects for a spin-2 
field (the graviton): radiation must be present only in the spatial tensor mode. Indeed, 
from the Einstein equations one can show that the tensor mode obeys the wave equation 

d 2 T + 2r]d T - V 2 + 2K) = 87rGa 2 S iiiT , (21) 

where rj = dloga/dlogr, V 2 = ^ViV j is the spatial Laplacian, K is the spatial curvature 
constant, and S^t is the proper transverse-traceless stress. 



The field equations fll9D-(j2"ID show that for a fluid with velocity v, W{ = ip x 0(v/c) and 
hij — tp x 0(v/c) 2 . These relations will be important when we compare with the Newtonian 
limit. 

It is straightforward, though algebraically tedious, to compute the components of the 
Weyl tensor for the metric given above. The result is: 

C° i0j = -~ A;(V> + 4>)~\ Wij + \{d 2 T + V 2 - 2K)h i3 , (22) 

C° ijk = {VkW^ - VjW ik ) + ^(V 2 + 2K){ lijWk - ^ kWj ) + (Vjkk - Vukj) , (23) 
C^u = l l ]mn l pn kl [C m 0p0 + (V 2 - 3K)h m p ] + l l ]mn l pq kl V p V n h m q , (24) 



- 9- 



where 7 U is used to raise the components of V, and hij and we have defined several auxiliary 
quantities: 

D i:i = V.Vj- - - 7ij V 2 , Wij = V (i w i} , 7 * ifcl = <f kljl - 5 l lljk ■ (25) 

All other components of the Weyl tensor follow from the symmetry relations 

C^ UK \ = C[^ u ][ K x\ = C K \^ LU . (26) 

Note that because we are retaining only first-order terms in the metric perturbations, the 
Weyl tensor components may be raised and lowered using the unperturbed metric. Thus, 
up to powers of the expansion factor, we may regard equations ([22]) -(|24|) as giving us the 
Weyl tensor components in locally flat coordinates. 

The electric and magnetic parts of the Weyl tensor follow from equations @. We first 
give them using the 4-velocity of comoving observers, = u% = (a _1 ,0) for the metric 



of equation (17). Note that because the Weyl tensor is already first-order in the metric 
perturbations, we need the 4-velocity only to zeroth-order. We obtain (cf. Bruni, Dunsby, 
& Ellis 1992, eqs. [113] and [115]) 

Etiiuc) = \ Ai(^ + 0) + \ Wij - \{d 2 T + V 2 - 2K)h l3 , (27) 

HiAuo) = ~\ VaHfl + V k e kl {i h 3)l . (28) 

These components are given in the comoving coordinates of equation (P^7|). To obtain the 
components in a locally flat coordinate system at rest relative to the comoving frame (to 
first order in the metric perturbations) one simply multiplies Eij and by a~ 2 . 

In the Newtonian limit, ip = <p is the Newtonian potential and is simply the 
gravitational tidal field (the traceless double gradient of the potential, with a factor of a~ 2 
required to convert the gradients from comoving to proper coordinates). The magnetic 
part of the Weyl tensor appears to have no Newtonian counterpart. It depends on 
gravitomagnetism and gravitational radiation, as do the corrections to the electric part. 

Next we must obtain the electric and magnetic parts of the Weyl tensor in the fluid 
frame moving with 4-velocity w M . Two steps are required for this computation. First we 
must evaluate E^ u and H^ u using u M = a -1 (7, 7^) rather than u£ in equation where 
v is the peculiar velocity and 7 = (1 — f 2 ) -1 / 2 . Then we must transform from comoving 
coordinates to locally flat coordinates in the fluid rest frame. The result is 

E'^u) = a- 2 A^(v)A^(v)u K u x C, KuX (29) 
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and similarly for H iy The Weyl tensor may be computed from Eij(u c ) and Hij(u c ) using 
equation (|7|). A prime is used to indicate that the components are given in the transformed 
frame moving with velocity v relative to the unprimed (comoving) frame, and A k i (v) 
is the Lorentz transformation corresponding to a boost v: A° = 7, A° i = A l = jv 1 , 
A* • = 5 l j + (7 — l)v l Vj/v 2 . We are allowed to use special relativity here because we are 
working in locally flat coordinates. [The factor a -2 converts E^iu) from comoving to 
locally flat coordinates as described in the previous paragraph.] 

The two-stage transformation described above gives the following results: 

% = E l3 + 2v k e kl (i H 3 y + 0(v/c) 2 , H' tJ = H tJ - 2v k e kl {l E 3)l + 0(v/c) 2 , (30) 

where, in each frame (primed and unprimed) the components are evaluated using locally 
Minkowski coordinates. (To reduce the clutter we have dropped the arguments u and u c 
from the primed and unprimed tensors, respectively.) Note that, although we began with 
a covariant definition of E^ and H^, here we are evaluating them in two different locally 
inertial frames such that in each frame the tensors are purely spacelike. Our results are 
reminiscent of the Lorentz transformation of the ordinary electric and magnetic fields. Ellis 
(1973) shows how the electric and magnetic fields may be defined as spacelike 4-vectors 
using the 4-velocity field to split the electromagnetic field strength tensor in a way similar 
to what we have done here for the Weyl tensor. 



4. Newtonian Evolution of the Weyl Tensor 

In the preceding sections we have derived the equations of motion for the electric 
and magnetic parts of the Weyl tensor using general relativity, and we have related these 
fields to the Newtonian gravitational potential and other metric perturbation fields of a 
perturbed Robertson- Walker spacetime. In this section we provide a Newtonian derivation 
of these results for cold dust. We will work in a perturbed Robertson- Walker spacetime, 
but a similar derivation can be carried through in a non-cosmological setting. 

We begin with the continuity and Poisson equations in comoving coordinates 
(Bertschinger 1992): 

85 

— +V- [{l + 6)v\ =0 , V 2 = 4vrGa 2 (p-p) . (31) 

The mass density is written p = p(r)(l + 5) and the peculiar velocity is v — dx/dr. 
Gradients are taken with respect to comoving coordinates. The cosmological Poisson 



- 11 - 



equation differs from its flat spacetime counterpart in that the mean background density is 
subtracted from the source. Equations (|3T| ) are valid to first order in v/c on scales much 
smaller than the Hubble distance cif -1 . Note that these equations are valid even after 
orbits intersect provided that 5 and v are taken to be the total density fluctuation and the 
average fluid velocity 

We define the Newtonian gravity vector in comoving coordinates, 

g = -V(j) , V • g = -AnGa 2 p 5 . (32) 



This differs from equation fl2"0|) in that we neglect d T w. For slowly moving sources we expect 
d T ~ v ■ V so that the neglected contribution is 0(v/c) 2 compared with the Newtonian 
gravitational field. 

For convenience we will decompose the velocity gradient into trace and traceless 
symmetric and antisymmetric parts as in equation (0): VjU,- = \05ij + <7y +tOij. Note that 
we use 9 = V ■ v for the divergence of the peculiar velocity using comoving coordinates; 
O = 9/a + 3H is the proper divergence of the proper velocity including Hubble expansion. 
We decompose the gravity gradient in a similar manner: 

- V, 9j = ViVj = i 6 {j V 2 + Ea . (33) 



This equation defines in the Newtonian case; note that here, as in equation (|2j 
the comoving components (Cartesian so that 7^ becomes 5ij), which differ by a factor a 2 (r) 
from the proper components. Equations ( P7[ ) and Q3"3] ) agree to 0(v/c). 

Returning now to our derivation, we substitute g into the continuity equation to obtain 
V- (j^f\ =4vrGa 3 V-/ , f = pv = f\\+h. (34) 

We have decomposed the mass current in the comoving frame into longitudinal and 
transverse parts obeying V X f\\ — and V ■ f± = 0. Integrating equation ([34]) yields 

f + ^^AnGa 2 ^. (35) 
or a 

The source involves only the longitudinal mass current because the Newtonian gravity 
vector is longitudinal. We rewrite equation ([35]) using the Lagrangian derivative following a 
fluid element, d/dr = d/dr + v ■ V: 

^l + ^g = y.y g + 4:1T Ca 2 h . (36) 
ar a 
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Before proceeding further we replace the transverse mass current by a transverse vector 
field H obeying the field equations 

V x H = -lQnGa 2 f ± , V • H = . (37) 

We recognize this as the gravitomagnetic field of equations ( |H?D and (POD- However, in the 
Newtonian case we imply no relation between H and spacetime metric perturbations; we 
simply regard H dynamical field related to f±. 

We further define a traceless tensor 

Hi, =~V j H i + 2v k e kl i V, 9l = Htj + e ijk A k , (38) 
which we have decomposed into a symmetric part 

H ij = -±V {i H j) -2v k e k \ i E j)l (39) 



and an antisymmetric part with dual 



Ai=^v { V 2 - Eijv* + \ ( V x H) . . (40) 

To first order in v/c, equation (^) agrees exactly with the magnetic part of the Weyl tensor 
in the fluid frame (eqs. [^8[] and |3~C)|| ). 

For reference we provide some useful identities following from the definitions of 
equations (|32|)-(|33D and (|37f)-(|40[). First are two differential identities for the divergence 
and curl of A: 

V t A { = -6\/ 2 (f>-^E lJ , (41) 

e kl i V k A l = -8nGa 2 p (l + 1 ^ - e ijk a jl E\ + E ij( J . (42) 

Next are two differential identities for Hij. The first is 

VjH'i + e ijk a jl E k l + ZE^J = -8nGa 2 pUj . (43) 

Aside from the factor a 2 present on the right-hand side because we are using comoving 
rather than proper coordinates, this result agrees exactly with equation (|15|). The last 
identity will be useful in deriving equation flip]): 



e kl 4 V fc ^- = -VjAi + - Vj-vt + - 9 Eij + S tj a kl E kl - Aa k (i E j)k - Au%E j]k . (44) 
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To derive equation flTBD we first substitute equation ( ^0|) into equation (|37|) and use 
equations fl3"2| ) and ( |36| ) to get 

~r + -g — A-kGcl pv + A . (45) 
get a 

Taking the gradient gives 

+ a - V, ^ + (V, u fc ) ( V, 9j ) = A7tGa 2 p V, ^ + V, A j . (46) 
The traceless symmetric part of this equation is 

—j^- + — Eij + -9 Eij + - (Tij V 2 + a k ^Ej) k — - 5^ o kl E k i — u h ^E^ h = 

- 47iGa 2 p on - V (i + ^ 5ij V ' k A k . (47) 



The trace and traceless antisymmetric parts of equation (46) give no further information 
because they simply reproduce the divergence and curl of A that are given above. 



Equation (|47D involves the symmetrized gradient of the antisymmetric part of Hij 
defined by equation (0), which can be replaced in favor of the antisymmetrized gradient of 
the symmetric part ify using equations flHD and fl44]). The necessary link is the identity 

VpA,-) - i % V fe A fc = - V fc e H + ~ ^ + H ^ . V 2 - 4<7* (< £ i)Jfe + | % a kl E k i . (48) 



When this is substituted into equation fl4"?|) we obtain 



^ + ^ £ y - V, e fe < ^ify, + 0£ y + u kl E k i - 3a k (i E j)h - iv h (i E j)k = -A-nGa 2 pa l3 . ( 19) 



If the variables are converted to proper coordinates, dt = adr, E^ — » a 2 Eij, etc., and one 
recalls that the comoving expansion scalar is 9 = aO — 3d/ a, one obtains equation fllTf ). 

We have succeeded in deriving the div-if and E equations from Newton's laws; 
what about the other two field equations for the Weyl tensor? Care is required because 
these equations involve terms 0(v/c) 2 , as we can see by the following argument. From 
their definitions, we see that the units of and E^ differ by one power of a velocity: 

~ vEij. From their respective field equations, one sees that v must be of order the 
matter peculiar velocity. Dimensional analysis of equations (|13|) and ([14]) shows that the 

terms must be divided by c 2 and that they are 0(v/c) 2 relative to the other terms. 
(The same analysis applied to eqs. |H| and [16j shows that all terms are of the same order.) 
In the strictly Newtonian limit one would drop terms 0(v/c) 2 ; for example, equations (|3T|) 
neglect relativistic corrections of this order. 
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Indeed, equations (13) and (14) are satisfied in the Newtonian limit, as one can see by 
evaluating the divergence and curl of Ey from equation fl3"3|): 

Vj&i = 2 - V, V 2 , V fc e kl {t E j}l = . (50) 

Using the Poisson equation for V 2 0, one obtains 

V 3 E\ = ^Ga 2 V tP . (51) 



Aside from the use of comoving coordinates, this equation and the second of equations ( |5"0"D 
agree with equations (|T3"| ) and (O) in the Newtonian limit. 

This derivation is not fully satisfactory because we have not accounted for all the terms 
in equations ( |T3|) and fll4|). Fortunately, with a little care we can derive the full equations. 

First, we must use the E^ and obtained using general relativity rather than 
Newtonian gravity, although we shall make the simplifying assumptions that gravitational 
radiation and the gravitational effect of shear stress, as well as terms that are explicitly 
0(v/c) 2 , can be neglected. Using equations ([27D, fl28|) , and (j30|), in comoving coordinates 
we obtain the results 

1 ~rA a , 1 Ti/. . i o-.. ,*l tr(°) tr - 1 V7. tr , o„.. M 



^tf = (V,V,- - - V 2 J + - W« + 2^6*' (i JT$ , Jfy = -- V (l ^) - 2^ e fe( ^ , (52) 

where superscript (0) indicates the quantity is to be evaluated setting t> fc = 0. We must 
include in equations ( [52D terms that are explicitly first-order in v^, even though we will 
evaluate the evolution equations in the fluid rest frame for comparison with equations ( |i~3"D 
and (|T4]), because we will need spatial derivatives of Eij and Hy. This fact shows that we 
can, however, safely ignore terms that are explicitly quadratic in v^. 

The divergence of E^ from equation fl52|) gives 

VjE\ = y Ga 2 V lP + ^V 2 w t + e ijk a jl H k t + 3H i:j uj j , (53) 

where we have retained only the terms that do not vanish when v k = 0. If we can show that 
V 2 Wi = d T (16nGa 2 f±i) vanishes in the fluid frame, then equation ( |53| ) implies equation 
(|T3|). We can show this using the following trick. In the fluid frame, f — f\\ + f± — 0, by 
definition. Also, the fluid frame is freely-falling, so that the fluid acceleration measured at 
the fluid element must vanish. Equation (j36D then yields f\\ =0, hence f± = and therefore 
V 2 Wi = 0. Since this result holds for all time (the fluid continues to remain at rest in this 
inertial frame), V 2 Wi = 0. Equation (0) then follows. 
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To derive the H equation we take the curl of E i3 - from equation (|52|), again discarding 
terms that do not vanish when Vk = 0: 

V k e kl (i E 3)l = \v {i H 3) - 6H i3 - 5 i3 a kl H kl + 3a k (i H j)h + u k {i H j)k . (54) 

In evaluating the time derivative term we must be careful to divide the comoving 
components H 3 used here by a to obtain the proper components. Then, using the fact that 
in the fluid frame v k = dv k /dr = 0, we obtain dH^/dr = d T Hi 3 = — (<9 T — d/a^V^-f/j). 
Equation (|54D then reduces to 



+ ~ Hij + ^ eH (*-%)* + 9Hij + 5ij a kl H k i — 3cr k ^Hj) k — Lu k ^Hj) k = . (55) 



Converting to proper coordinates and using the relation 9 = aO — 3d/ a, we recover equation 
(0)- 

It is interesting to note that in the Newtonian derivation of the iJ-dot equation (14) 
- but not the £"-dot equation (16) — we had to assume vanishing acceleration in the 
fluid frame. The reason is clear from the relativistic equations (9) and (11) retaining the 
4-acceleration. By dimensional analysis, the 4-acceleration term in equation (11) must 
be divided by c 2 compared with the term in equation (9). Thus, acceleration effects are 
unimportant in the Newtonian tidal evolution equation but are important in the evolution 
of H i:j . 



5. Conclusions 



We have shown the magnetic part of the Weyl tensor Hij does not necessarily vanish 
in the Newtonian limit by deriving an explicit expression for it — equation (39). As a 
consequence, the Lagrangian evolution of the tidal field (the electric part of the Weyl 
tensor, Ey) is not purely local: it involves the gradient of (eq. [16|). This term is 



of the same order in powers of v/c as the other terms affecting the tidal evolution. No 
ambiguity remains because we have obtained identical results using both general relativity 
and Newtonian gravity. 

It may be surprising that magnetic-like effects are significant for the evolution of 
the tidal field. After all, the Newtonian limit neglects the gravitomagnetic force mv x H 
relative to the gravitoelectric force mg. However, changes in the tidal field are dependent 
on the motion of the fluid. As in the case of electromagnetism, electric fields in one frame 
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transform into magnetic fields in another frame. Moreover, even if we choose a frame so 
that Hij = (not necessarily the fluid frame), this condition is insufficient to restore locality 
because the evolution of Eij depends on the gradient of H^. Given the complexity of the 
equations, at this time we cannot provide more physical insight into our results. However, 
it is clear that they are implied by mass conservation and Newtonian gravity rather than 
some subtle relativistic effects. 

Although we have shown that the tidal evolution is, in general, nonlocal, the Lagrangian 
fluid method still can be applied to study the nonlinear evolution of self-gravitating mass. 
One will simply have to evolve both and and compute their gradients. Since this 
requires knowing the positions of mass elements, the trajectories will have to be integrated 
simultaneously. Although this greatly complicates matters, there may be advantages in 
being able to compute the density and velocity gradients of individual mass elements 
treated as a fluid rather than by summing over particles weighted with a smoothing kernel. 
After trajectories intersect, however, the Lagrangian fluid method becomes substantially 
more complicated. 

Another possibility is that an alternative approximation may prove useful in restoring 
locality. Since we have obtained an explicit Newtonian expression for H^, one can test 
alternative local approximations for its gradient. Of course, we still have not shown how 
bad an approximation it is to neglect altogether. This can only be done by comparing 
with other approximations or N-body simulations. This work is currently underway. 

Our results do not resolve the question of whether nonlinear gravitational collapse 
favors the formation of prolate filaments as opposed to oblate pancakes. This question 
hinges critically on the behavior of during gravitational collapse. We leave this subject 
for future work. 
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